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ABSTRACT 


In  a  previous  paper  (“A  Multivariate  Exponential  Distribution," 
Boeing  Scientific  Research  laboratories  Document  Dl-82-0505)  the  authors 
have  derived  a  multivariate  exponential  distribution  from  points  of 
view  designed  to  indicate  the  applicability  of  the  distribution.  Two 
of  these  derivations  are  based  on  "shock  models"  and  one  is  based  on 
the  requirement  that  residual  life  is  independent  of  age. 

The  practical  importance  of  the  univariate  exponential  distribution 
is  partially  due  to  the  fact  that  it  governs  waiting  times  in  a  Poisson 
process.  In  this  paper,  the  distribution  of  joint  waiting  times  in  a 
bivariate  Poisson  process  is  investigated.  There  are  several  ways  to 
define  "joint  waiting  time."  Some  of  these  Jead  to  the  bivariate 
exponential  distribution  previously  obtained  by  the  authors,  but  others 
lead  to  a  generalization  of  it.  This  generalized  bivariate  exponential 
distribution  is  also  derived  from  shock  models.  The  moment  generating 
function  and  other  properties  of  the  distribution  are  investigated. 


1. 


Introduction . 


It  is  common  in  applications  to  encounter  random  vectors  with 
elements  that  are  clearly  dependent.  In  such  cases ,  the  marginal 
distributions  are  often  known,  but  these  alone  do  not  determine  the 
joint  distribution;  to  the  contrary,  it  is  well  known  that  the  Joint 
distribution  can  take  many  different  forms.  However,  it  may  be  that 
the  joint  distribution  can  be  obtained  through  a  study  of  the 
mechanism  causing  dependence.  In  the  case  of  exponential  marginals, 
we  have  in  a  previous  paper  (Marshall  and  Olkin,  1966)  presented  such 
derivations  of  the  bivariate  distribution  given  by 

(1.1)  P{X  >  x,  Y  >  y]  =  F(x,y)  -  for  x%y  >  0> 

If  X  and  Y  are  life  lengths  of  devices  subjected  to  shocks,  this 
distribution  arises  when  the  occurrence  of  shocks  is  governed  by  one 
or  more  Poisson  processes.  In  addition,  (1.1)  is  the  unique  bivariate 
distribution  (with  exponential  marginals)  which  satisfies  the  require¬ 
ment  that  residual  life  given  survival  to  a  common  age  t  has  a  dis¬ 
tribution  independent  of  t. 

The  viewpoint  of  this  paper  is  suggested  by  the  fact  that  the 
univariate  exponential  distribution  governs  waiting  times  in  a  Poisson 
process.  More  precisely,  if  { Z(t),  t  i  0}  is  a  Poisson  process 
with  parameter  \  and  s  is  a  fixed  point,  then  the  waiting  time  X 
from  8  to  the  next  event  (jump)  of  the  process  has  an  exponential 
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distribution  with  parameter  X.  Likewise,  if  s  is  replaced  by  a 

random  variable  5,  the  occurrence  time  of  the  i ^  event,  then  the 

s  t 

waiting  time  X  from  S  to  the  t+1  event  is  again  exponentially 
distributed  with  parameter  X.  More  generally,  S  can  be  any  stop¬ 
ping  time,  i.e.  ,  any  non-negative  random  variable  for  which  one  can 
check  whether  S  <  s  by  observing  X(t)  only  for  t  <  s.  To 
avoid  ambiguities  and  problems  of  separability,  we  assume  throughout 
this  paper  that  sample  functions  of  Poisson  processes  are  right  con¬ 
tinuous  with  probability  one. 

Following  Dwass  and  Teicher  (1957),  we  define  (in§2.1)  a  bivariate 
Poisson  process  { Z(t)  *  (Z; (t) ,Z2(t))  ,  t  1  0}  where  Z\(t)  and 
Z2(t)  are  correlated  Poisson  processes.  There  are  several  ways  to 
define  a  "joint  waiting  time"  for  this  process: 

( i )  We  may  choose  a  fixed  time  s,  and  consider  the 

waiting  time  X  to  the  next  event  in  the  Z\  process 

together  with  the  time  Y  to  the  next  event  in  the 

* 

Z2  process. 

(li)  In  place  of  a  fixed  time  s  as  in  (£) ,  we  may  gin 

waiting  at  a  random  time  S  that  is  a  stopping  time. 

(Hi)  We  may  consider  the  waiting  time  X  from  a  fixed  point 

s  to  ne  next  event  in  the  Z\  process  together  with 

the  waiting  time  Y  from  the  fixed  point  s+6  to 
the  next  event  in  the  Z2  process.  Alternatively,  it 
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may  be  that  replaced  by  a  stopping  time. 

( iv )  F  nation  is  obtained  if  in  (Hi) ,  both 

8  and  s+6  are  stopping  times. 

It  is  not  difficult  to  show  that  in  cases  (i)  and  ( ii )  ,  the  joint 
distribution  of  X  and  Y  is  given  by  (1.1).  On  the  other  hand, 

(Hi)  leads  uo  the  generalization 

(1.2)  P{X  >  x,  Y  >  y}  =  F(x,y;6)  - 

■  e-A1ar-A2!y“A12max[x,i/-hnin(a:,6) )  g  >  Q 
e-\lX-\2y-*i2max[x,y+min(x,-6)]'  g  <  o,  xty  >  0. 

The  parameter  6  might  be  called  a  "shift"  parameter,  though  it  is 
not  simply  a  location  parameter.  Of  course  when  6  ■  0,  (1.2) 
reduces  to  (1.1). 

Case  (iv)  above  results  in  a  still  more  general  family  of  dis¬ 
tributions.  These  can  be  obtained  from  (1.2)  by  first  conditioning 

on  6,  so  that  the  result  is  a  mixture  of  such  distributions: 

00 

(1.3)  P{X  >  x,  Y  >  y]  -  /  F(x,y;6)dG(6). 

—00 

Because  of  this  representation,  we  confine  ourselves  in  this  paper 
almost  exclusively  to  a  discussion  of  the  distribution  given  by  (1.2). 

In  addition  to  governing  the  waiting  times  in  the  bivariate 
Poisson  process,  (1.2)  can  be  derived  from  shock  models  (§2.3). 

These  derivations  shed  light  on  the  applicability  of  (1.2),  particularly 
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in  reliability  problems. 

For  convenience,  we  say  that  X  and  Y  are  bivariate  exponential, 
BVE(Ai ,A2 »^12 »fi) *  if  (1.2)  holds,  and  refer  to  the  distribution  of 
(1.2)  as  the  BVE(Ai , A2 , Aj2;6)  distribution.  When  6  is  the  only 
significant  parameter,  it  is  convenient  to  abbreviate  BVE(Ai,A2, Ai2;6) 
by  writing  BVE(6).  We  refer  to  F(x,y;6)  =  P{X  >  xt  y  >  y]  as  the 
survival  probability,  this  has  a  much  simpler  form  than  the  distribution 
function  F(x,y\&)  =  P{X  f  x,  Y  <  y\. 

The  survival  probability  F(x,y;6)  is  clearly  decreasing  in 
6  2  0  and  increasing  in  6*0.  Since  F(x,z/;6)  has  marginal  survival 
probabilities 

Pl*  >  x)  -  nx.0-,6)  -  e‘(>1+>12)* 
p{y  >  y}  °  P(0,^;6)  *  e  ^2+^12 )i/ 

that  are  independent  of  6 , 

F(x,y\6  j)  -  F(x,y  ;62)  -  F(x,z/;  6j)  -  F(x,z/;62). 

Thus,  F(x,y\ 6)  is  also  decreasing  (increasing)  in  6^0  (6*0). 

If  F(x,z/;6)  is  the  survival  probability  of  a  pair  of  items,  it 
follows  that  the  probability  both  items  survive  and  the  probability 
neither  item  survives  both  reach  a  maximum  when  6*0  (the  case  of 
(1.1)),  and  a  minimum  when  6  ■  ±°°  (the  case  of  independence). 

Various  properties  of  the  distribution  F(x,y’,6)  are  discussed 
in  §3,  generalizations  are  obtained  in  §4,  and  a  brief  discussion  of 
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the  multivariate  case  is  given  in  §5. 

2 .  Derivations. 

2.1.  The  bivariate  Poisson  process. 

The  univariate  Poisson  distribution  can,  of  course,  be  obtained 
as  a  limit  from  the  binomial  distribution,  which  in  turn  arises  via 
convolutions  from  the  Bernoulli  distribution.  This  derivation  can 
be  repeated  in  the  bivariate  case  to  obtain  a  bivariate  Poisson  dis¬ 
tribution,  because  there  is  a  unique  bivariate  distribution  with 
Bernoulli  marginals.  That  distribution  places  mass  only  at  (0,0), 
(0,1),  (1,0)  and  (1,1).  By  convolving  this  distribution  and  taking 
limits,  Teicher  (1954)  obtained  the  bivariate  Poisson  distribution 


(2.1)  P{Z1  -  x,  Z2  -  y } 


.-»<P00«Pirt> oi)a%*’l'>  (Xpoo)°(Xpol)*~tt(Xpio)ir° 

a  a!  (x-a)  !  (y-a)  ! 


Earlier,  this  distribution  was  obtained  by  McKendrick  (1926),  who 
obtained  (2.1)  as  the  solution  to  a  difference-differential  equation. 
Starting  with  the  bivariate  Bernoulli  distribution,  Campbell  (1934)  obtained 
(2.1)  via  generating  functions. 


Dwass  and  Teicher  (1957)  have  shown  that  this  distribution  is  the 
only  bivariate  distribution  with  Poisson  marginals  that  is  infinitely 
divisible.  Now  if  { Z\(t)  t  t  1  0}  and  {Z2(£),  t  1  0}  are  Poisson 
processes  and  if  {(Zj(t),  Z2(t)),  t  >  0}  has  stationary  and  independent 
increments,  then  the  increments  must  have  an  infinitely  divisible  dis¬ 
tribution,  hence  the  distribution  given  by  (2.1).  Because  of  this 
uniqueness,  we  call  the  process  the  bivariate  Poisson  process. 


Dwass  and  Telcher  have  derived  the  bivariate  Poisson  process 
{(Zi(£),  ^2 ( t) ) »  t  1  0}  directly  from  a  one-dimensional  Poisson  process 
t  1  0}  as  follows:  If  Z*(fc)  jumps  at  time  t,  then 
Z(t)  ■  Z2<t))  jumps  at  t  from  Z(t-)  to  Z(t-)  +  (t,j)  where 

-  (0,0)}  -  pn,  P{(t,j)  -  (0,1)}  -  p10,  P{(i,j)  -  (1,0)}  -  p01, 
and  P{(i,j)  ■  (1,1)}  ■  Poo*  The  increments  (t,j)  at  each  jump  are 
independent  and  {Z*(t),  t  >  0}  has  parameter  X,  so  that 


PfZjU)  -  x,Z2(t)  -  y} 


mln(x,y) 

£  > 


L  : 

•x+y- a 


~Xt(Xt)k  klP<>0  p01  _ 

k !  a!  (x-a)  !  {y-a) !  (k-x-y+a) ! 


x-a  y-a  k-x-y+a 

p  10  p„ 


■  e 


-Xt(pnn4t,,n-tP„)'ll^)(XtP00)a(XtP0I)I'a(>*PlO) 

~0  a!  (x-a)  1  (y-a) ! 


y-a 


Directly  from  the  model,  it  is  easily  seen  that 

-Xt(poo+P01>  _ 

P{zx(t)  -  x}  -  - jj -  [At(p004poi)]  . 

In  their  derivation,  Dwass  and  Teicher  take  pn  ■  0,  but  this 
results  in  no  loss  of  generality. 

One  can  also  obtain  the  bivariate  Poisson  process  from  three  inde¬ 
pendent  one-dimensional  Poisson  processes  { Zf(t),  t  >  0},  |Z*(t),  t  >  0} 
and  {2*2(t),  t  >  0}  with  respective  parameters  Xj,  X2  and  Xi2.  In 
fact,  { (Z*(t)  +  ZJ2<t),  Zj(t)  +  2*2(t)),  t  >  0}  is  a  bivariate  Poisson 
process.  This  observation  is  utilized  in  12.3. 

2.2.  The  joint  waiting  time. 


Choose  s  j  and  82*  and  consider  the  time  X  from  s1  to 
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i 

I 

!  the  next  event  in  the  process  \Z±(t) ,  t  -  0}  together  with  the  time 

Y  from  8 2  to  the  next  event  in  the  process  { Z2 (t)»  t  >  0}.  To 
I  find  the  joint  survival  probability  P{X  >  x,  Y  >  yj ,  it  is  convenient 

to  assume  first  that  6  ■  82  -  -  0.  Let  -  8\  m  x  and 

t2  ”  s2  *  y-  The  three  cases  (t)  x  <  6,  (ti)  {  <  x  i  1/  +  {, 

1 

and  (tii)  zy  +  6  i  a;  (indicated  in  Figure  1)  must  be  treated  separately. 


y 


Case  (t):  x  <  6  (fll  i  tj  i  «2  i  t2)> 


p{*  >  x,  y  >  zy}  - 

•  • 

.  AX(t1-81)f  .  0  v[A(t2-S2)]J 

*  <4l-8l) - IT— - 71 - (Pll+p01)J 


j! 


»  e-*(Poo-*Poi)*-x(Poo+Pio)*/ 
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Case  (ii):  S  <  x  <  y  +  6  (si  1  e2  f  t i  i  t2) . 

P{X  >  x,  y  >  y}  -  2,  e  X(*2  8>) - Tl - (pu-^10)  • 

i.i.k 

.  e“APOO<S"^P01aj"^(POO'^10)i/t 

Case  (ttt):  z/  +  6  <  x  (a x  <  s2  -  *2  -  *i)  • 

r  i  V  -A(so-s i )  ^82”8 ^  J  * 

P{x  >  x,  y  >  I/}  -  2  e  A(82  8l) - -Tj - 

i.*.* 


'(Pll+Pio)' 


.  -A(t2-a9)tX(t2“82)1  .7  -A(ti-t9)[X(tl“t2)i 


J! 


Pile 


Xi  <Pir*Pio) 


“ACpoo^Poi^^Pioi/ 


Letting  Aj  ■  Apoi*  A2  *  Apio*  80(1  A12  “  Ap 0o»  we  obtain  from  these 
three  cases  the  survival  probability  for  6-0.  If  6  i  0,  the 
result  is  then  obtained  by  interchanging  Aj  and  A2,  x  and  yt 
6  and  -6.  Together,  this  yields 


(2.2)  P{X  >  x,  X  >  y] 


r  “Aix-A2&-A12(arty) 

e  1 

{  0-A1x-A2y-A12(x-6) 

e  > 

<  -Ajx-A^-A^ 

C  » 

-AiX-A2i/-A12x 

C  • 

-A  iX-A2j/-A12  (z/+6) 

C  > 

-A!X-A2i/-A12(aH-t/) 
v  e  » 


6  i  -y 

-y  -  &  t  min(x-z/,0) 

min(x-i/,0)  <  6  <  0 

0  5  5  i  max(x-z/,0) 

max(x-i/,0)  <  6  <  x 

6  -  x, 


which  can  be  rewritten  as  (1.2). 
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i  2.3.  A  shock  model  derivation. 

We  have  previously  us  both  fatal  and  non-fatal  shock  models 
|  to  derive  the  BVE(Ai,A2,Ai2»0)  (Marshall  &  Olkin,  1966).  Such  models 

can  be  modified  to  yield  the  BVE(Ai,A2,Aj2*6)  for  any  6;  only  the 
fatal  shock  case  is  considered  here. 

i 

i 

< 

Suppose  that  a  device  is  placed  in  an  environment  where  it  is 

i 

,  subject  to  fatal  shocks  from  two  sources,  governed  by  Poisson  processes 

^  Zj(t),  Z^2 (£)  with  parameters  A},  A 12  •  After  a  time  5,  a  second 

,  device  is  placed  in  a  similar  environment  where  fatal  shocks  are 

governed  by  the  processes  Z2(t),  Z^C t)  with  parameters  A2,  A12. 
Suppose,  in  addition,  that  the  three  processes  are  independent. 

Let  X  ( Y )  be  the  age,  or  length  of  service,  of  the  first  (second) 
device  at  the  time  of  death.  To  obtain  the  joint  survival  probability 
p{x  >  x,  Y  >  z/},  we  again  consider  the  cases  represented  by  the  three 
regions  in  Figure  1. 

Case  ( i ) :  x  *  6. 

i 

P{X  >  X.  y  >  {,}  -  .-Wl+»12>*-a2+»12>» 

/ 

Case  ( ii ) :  6  <  x  1  y  +  6. 

p{x  »  X,  y  >  y)  -  e-»i*-»2i/-»i2(irM) 

Case  (Hi) :  y  +  6  i  x. 

P{X  >  x,  i  >  y]  ■  e-»l*-»2ir»l2*. 


Thus  (XJ)  has  the  joint  distribution  given  by  (1.2). 
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3.  Properties  of  the  distribution  F(xty; 5). 

3.1.  The  distribution  function. 

The  BVE  distribution  (6-0)  has  both  an  absolutely  continuous 
and  a  singular  part,  with  positive  mass  on  the  line  x  *  y.  From 
the  shock  model  of  §2.3,  we  see  that  if  (X, Y)  is  BVE(6),  then  in 
case  an  event  in  the  process  Z\2(t)  causes  both  devices  to  fail, 

X  ■  X  +  6.  Thus  the  BVE(6)  distribution  has  a  singular  part,  with 
positive  mass  on  the  line  x  •  y  +  6. 


Theorem  3.1.  If  F(xty )  is  BVE(Xj  t\i  ,Xi2»$) »  <5-0,  and 
A«Al+A2+*12»  then 

F(x,y)  -  a F  (x$y)  +  a  F  (x,y) , 
a  o 


where 


oFa(xty) 


a 


1 


“(Ai+Ai2)6 

X  e 


a 


1  -  a, 


e-Aia>A2J/-Ai2inax[£,z/+min(x,6)  ] 


^12  A2<5-Amax(a:,z/+6) 

X  6 


is  absolutely  continuous,  and 

Fs(x,y)  -  e-^ [max(x,zy+6)-5 ] 

is  a  singular  distribution. 


Proof. 


Let  ®i  "  Xj  +  X j2*  ®2  "  X2  X 12 *  Since 


F(u,v) 


e 

e 

e 


-6jir  62V 

> 

—  A @2  A 
-SjU-X 2V 


*  <  6 

6  <  u  <  &  +  V 

6  +  v  £  u. 


-li¬ 


lt  follows  by  differentiation  that 


a fa(u,v)  - 


0ie2e_ei“"O2t’, 

u 

<  6 

6 

<  u  <  6  +  t; 

e1A2e-eiu-i2V, 

6 

+  v  <  u. 

By  computing  oF^x^y) 
the  condition  F  (0,0)  -  1. 
aF(xty)  -  F(xty)  -  oF(xty) 

o  CL 


l 


du  fdv  a /  (u,y),  we  can  obtain  a 


Jy 


a 7  ( xty )  is  obtained  by  subtraction: 
8 


from 


In  order  to  carry  out  this  program,  we  begin  by  computing  three 
integrals  which  can  be  combined  to  give  a Fa(x,y) .  These  integrals 
again  correspond  to  the  regions  shown  in  Figure  1. 


Case  (t) :  x  1  6. 


Let  A(xty)  be  the  integral  of  a fa  over  the  region 
{(u,u):  x  <  u  <  6,  V  >  y).  Then 

A(xty)  ■  Jdu  Jdv  a fa(utv)  ■  e  02^[e  0ia:-e  01 
x  y 


Case  ( ii ) :  6  <  x  <  y  +  6. 

Let  B(x,y)  be  the  integral  of  af  over  the  region 
{(u,y):  x  <  u  i  v  +  6,  v  ±y}>  Then 

B(x.y)  -  fdv  fiu  af  lu.v)  -  . 

y  x 
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When  x  -  y  +  <$, 

SQ/+6,#)  =~e"6l6"^; 

when  x  «  6, 

5(6,1/)  -  e'0*6"^  _  lL  -0i 6-ty 

A 

Case  (Hi):  x  >  y  +  &, 

Ut  C(X,!/)  be  the  ‘nigral  of  ofa  over  the  region 
((n,v):  u>x,  y<v<u-  6}.  Then 

,°°  w-6  0 

m  du  I  dv  af  (u,v)  -  e"A25-0l*  .  JL  A26-A* 

J  J  Cl  A  e  * 

*  */ 

When  x  «  y  +  6, 

<?<i/+6.!/)  -  ~  e’^'0!6. 

By  combining  the  various  integrals,  we  obtain  a F  • 

a ' 

Case  Lid  :  x  <  <$. 

+  B(6,y)  +  C(y+6,y )  «  e~el*-02Z/  _  lii  g-Sjfi-ty 

A 

g3se  (H) :  6  1  x  <  zy  +  <s . 

aFa(x»y)  -  5(x,i/)  +  C(y+6,y )  =  e~A l^-A i26-02z/  _  ,i~616-Azy 

A  “ 

x  >  1/  +  6. 

aFa(x*y )  *  5(x,x-<5)  +  C(x,^)  =  e~A2#-0i*  _  -Ax+A26 

A  e 


-13- 


Thls  completes  the  derivation  of  F  . 

a 


Note  that 

a  -  aF  (0,0)  -  1  -  e“(Xl+X12> <5 # 

C4  A 

It  is  now  easy  to  obtain 

(l-a)Fg(x,2/)  «  F(x,y )  -  aFa(x,y)  -  •—  e~X  max(x»y+i>)  +  x26s 

so  that 


F8(*>y) 


s-X[max(x,y+6)-6] 


3.2.  The  moment  generating  function. 

As  in  the  case  of  the  BVE  distribution  (6*0),  direct 

00  00 

computation  of  an  integral  of  the  fora  /  /  G(x,y)dF(x,y)  requires 

Jo  J  0 

separate  consideration  of  the  absolutely  continuous  and  singular  parts 
of  F.  However,  this  can  be  avoided  by  choosing  a  kernel  G  such 
that  G(0,y)  e  0  =  G(x, 0)  where  G  is  of  bounded  variation  on  finite 
intervals.  It  then  follows  from  integration  by  parts  (Young,  1917) 
that 

00  00  00  00 
J  f  G(x,y)dF(x,y)  mJJ  F(x,y)dG(x,y) . 


0  0 


0  0 


The  kernel  G(x,y)  ■  (1-e  )(l-e  s)  has  the  required  property  and 

in  addition,  the  integral  j J~GdF  is  a  moment  generating  function.  Thus 
we  compute  (again,  take  5  >  0) 

00  oo 

00  00 

4>(e,£;6)  =  £  J^(l-e~8X)  (l-e~ty)dF(x,y ;  5)  ■  st^J  F(x,y)e  8X~ty  dxdy 


by  breaking  up  the  integral  into  three  parts. 
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6  "  6  » 

fdxfa,  Hx.yU-'*-**  -fdxjfl  *-*l*-^-*l2  (^y)-ax-ty 


l  _  e“(A1+A12+«)4 
(  A  i+A  1 2+® )  (A2+A}2+t) 

oo  1/+6  00  w+6  .  .  _  t 

Jdyf*  dxnx,y)S,x-t»-f  dyf  dr  e-M^2ir*12<«H*>-**-tj, 

e“(Ai+Ai2+s)6 

(  A2+A i2+£)  (  A 1+A2+A  j  2+8+£) 

f  dyf  dxHx.yW™-^  «  f  dyf  dx  ^l^Zh^n^sx-ty 

0  J  y+&  0  #+<$ 

e“(A 1+A 12+s) 6 
( A j+A  i2‘*’s )  ( A 1+A2+A  i2+8+£) 


Summing  these  three  integrals  yields,  for  6  £  0, 


u.i)  iisjiiil. 

S  V 

Interchange  of  e 
for  6  <  0. 


(Ai+A^2+s)  ( A2+A 12+^)  (_ 
and  t,  A 1  and  A2,  6 


_  e~(  A  i+A  12+s)  6 

A+s+t  J* 

and  -6  yields  (j )(s,t)/st 


To  obtain  the  moments  of  the  BVE(6),  we  note  first  that  since  the 
marginal  distributions  do  not  depend  on  6, 


Ex  ■  life  ^  '  *  TIT^F 

E1!  "  *2+*12*  ^  7  ’  <*2+»12>2  ‘ 
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We  also  have,  for  6*0, 

Aioe“^Xl+Xl2*6 

EXY  ■  lim  iiSaA 1  .  _ i _ +  - 

8  ,£*K  st  (  A  j+A  ^2)  (  A2+A  1 2)  (A1+A12)  (  A2+A12)  A  * 

X12e-(>1+>12>4 

(3.2)  Cotttf.I)  -  X(X1+A12)(A2+X12) 

and  the  correlation 

(3.3)  Corr(XtY)  =  p (XtY)  -  ^  e‘(Xl+Xl2)6. 

Note  that  0  1  p(XtY)  f  1,  and  that  X  and  Y  are  independent 
when  p(Af,y)  ■  0. 


3.3.  Representation  In  terms  of  Independent  random  variables. 

Theorem  3.2.  For  6  >  0,  ( X,Y )  is  BVE( Ax , A2 , A12 ;5)  if  and  only 

if  there  exist  independent  exponential  random  variables  U i,i/2,t/3  and 
i/4  with  respective  parameters  Alt  A2  » A *^12  such  that 


X  «  min(i/1,i/3). 


y 


min(i/2,£/3-6), 

min(i/2,i/4). 


if  l/3  >  5. 

if  i/3  <  6. 


This  theorem  can  be  obtained  directly  from  the  shock  model  of  §2.3,  or 
it  can  be  verified  formally  from  the  relation 

p|*  >  x,  y  >  y]  -  pji/j  >  x,  U2  >  z/}[Pii/3  >  Xt  U3  >  y  +  6|i/3  >  5}p(i/3  >  6} 

+  P{i/4  >  i/}p{i/3  >  x\U3  <  5}p{i/3  <  6}]. 

This  characterization  may  be  of  some  interest  for  6  >  0,  even 


though  the  simplicity  of  the  case  6*0  is  lost. 
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In  case  6  ■  0  the  representation  X  »  min(£/i  ,^3) ,  Y  ■  min(£/2»^3) 
immediately  yields  the  f-^ct  that  min(Y,Y)  -  min^i,^*^)  is  exponentially 
distributed.  However,  if  (Y,Y)  is  BVE(6)  and  6  >  0,  then  n.  ( XtY ) 
is  not  exponentially  distributed,  but 

P(mln(jr,r)  >  u)  -  e-to-»12"ln<w.«>, 

where  A«Ai+A2+^12*  This  distribution  is  piecewise  exponential, 
and  has  a  decreasing  hazard  rate. 

4.  More  general  bivariate  exponential  distributions. 

In  §2.2  we  considered  the  joint  distribution  F(x,i/;6)  of 
waiting  times  X  (Y)  from  Sj  (s2)  to  the  next  event  in  the  Z\  (Z2) 
process.  As  mentioned  in  the  introduction,  s  1  and  e 2  can  be 
replaced  by  random  variables  (stopping  times)  S\  and  S2>  so  l°n8 
as  52  -  *  6  is  retained  as  a  fixed  number.  That  the  resulting 

waiting  time  distribution  is  unchanged  can  be  seen  in  a  variety  of 
ways,  e.g.  by  observing  that  F(xty\&)  depends  on  8^  and  s2  only 
through  6. 

If,  on  the  other  hand,  6  is  replaced  by  a  random  variable  A, 
then  the  joint  waiting  time  distribution  is  a  mixture  over  6  of  dis¬ 
tributions  F(x,y;(>)  (mixed  according  to  the  distribution  of  A).  Of 
course,  such  a  distribution  has  exponential  marginals,  because  the 
marginals  of  F(x,y;6)  are  independent  of  6.  If  A  -  S2  **  has 
the  right  continuous  distribution  G  and  X  (Y)  is  the  waiting  time 
to  the  next  event  after  Si  (52)  in  the  Z\  (Z2)  process,  then 
we  obtain  from  (2.2)  that 
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<4.i)  ■  f  <&<«>  - 

e-\ix-\2y  e-Aix-X2 y 

-  +  e->>^/”ln(X'i/,0)+eX>25IiC(4) 

-y 

+  e  *12^[G(0)  -  G(min(:c-z/,0))  ]  +  e  *12X[(?(max(:e-i/,0))  -  (7(0)] 

+  e~Xn6dG(&)  +  e^^Gix). 

max(x-y ,0) 

An  Interesting  special  case  of  (4.1)  is  obtained  with  5}  (S2) 
the  time  of  first  event  in  the  Zi  ( Z2 )  process.  These  are  natural 
times  to  initiate  waiting  periods,  and  one  might  hope  that  for  this 
case  (4.1)  would  take  a  relatively  simple  form.  Unfortunately,  this 
is  not  the  case.  To  see  this,  we  note  first  that  ( S\,S2 )  have  the 
joint  distribution  F(xty; 0).  Thus  A  ■  S2  -  Si  has  the  distribution 
G  which  with  0!  -  Xj  +  X12,  02  ■  X2  +  X12,  X  «  \\  +  X2  +  X12  is  given  by 

P{S2  -  >  6}  -  f  Xie~QzS~Xaidsi  «  Xie-e2<5/X,  6  >  0, 

Jo 

P{S2-Si  <  6}  -  P{Si  -  S2  >  -6}  -  X2e0l6/X,  6  <  0, 

P{52  —  5^  ■  o}  “  X12/X. 


It  follows  from  this  and  from  (4.1)  that  if  X  (I)  is  the  waiting  time 
between  the  first  and  second  events  in  the  Z\  (Z2)  process,  then 
for  x  <  y. 
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P{ X  >  x.  Y  >  u)  e 

02 ( X+X i 2 )  ^  X ! X ! 2  _(09+Xl9)x 

e-XjX-X2iy  X 

62+Xi2  82+A]2 

XzXn  -81  (y-x)  _  -(61+X12)x 

6i+X12  L 


5 .  Multivariate  exponential  distributions. 

The  extension  of  the  BVE(6)  distribution  to  higher  dimensions 
is  best  understood  by  reference  to  the  shock  model  derivation  of 
§2.3.  Suppose  that  device  i  is  placed  in  service  at  time  6^,.  For 
each  non-empty  subset  of  the  variables,  there  is  a  Poisson  process 
governing  shocks  to  members  of  the  subset. 


Let  S  =  {s  *  )  :  s.  =  0  or  1,  but  s  j  (0,...,0)}. 

For  each  8  e  S,  let  Zg  be  a  Poisson  process  with  parameter  Ag , 
and  suppose  that  an  event  in  this  process  corresponds  to  a  fatal  shock 
to  component  i  if  and  only  if  s.  *  1.  Suppose  that  the  processes 

*\s 

Z  are  mutually  independent.  Let  X.  represent  the  life  of  the 

S  'V 

i ^  component,  and  let  I.  be  the  interval  [5.,x.  +  6.].  Then  if 

"V  Z-  Z  Z 

n  is  Lebesgue  measure, 


(5.1) 


xi 


n 


>  x  } 
n‘ 


-£V(U<i:a.. 

sc  S  % 


i*1 


One  can  also  write  (5.1)  expli  tly,  giving  its  form  for  various 
regions  defined  by  inequalities  on  the  6-  and  x..  Assuming 

Z  Z 

6j  =  0  <  62  <  ■  *  <  6  ,  there  are  3*5 *7  •••  • (2n— 1)  such  regions, 


-I9¬ 


60  this  is  prohibitive  even  in  three  dimensions.  (5.1)  can  also  be 
written  using  maxima  and  minima,  but  this  also  is  unattractive  even  in 
three  dimensions.  On  the  other  hand  (5.1)  in  its  present  form  is 
both  compact  and  easily  evaluated  for  any  given  6.  and  x.. 

If 


-21- 


REFERENCES 

[1]  Campbell,  J.  T.  (1934).  The  Poisson  correlation  function. 

Proa.  Edinburgh  Math.  Soa.t  Series  2,^18-26. 

[2]  Dwass,  Meyer  and  Teicher,  Henry  (1957).  On  infinitely  divisible 
random  vectors.  Ann.  Math.  Statist.  28  461-470. 

[3]  Marshall,  Albert  W.  and  Olkin,  Ingram  (1966).  A  multivariate 
exponential  distribution.  Boeing  Scientific  Research  Laboratories 
Document  Dl-82-0505.  To  appear  in  J.  Amer.  Statist.  Soo. 

[4]  McKendrick,  A.  G.  (1926).  Applications  of  mathematics  to  medical 
problems.  Proa.  Edinburgh  Math.  Soa.  44  98-130. 

[5]  Teicher,  Henry  (1954).  On  the  multivariate  Poisson  distribution. 
Skand.  Aktuarietidskr.  37  1-9. 

[6]  Young,  W.  H.  (1917).  On  multiple  integration  by  parts  and  the 
second  theorem  of  the  mean.  Proa.  London  Math.  Soc.t  Series  2, 


16  237-293. 


